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OO . ABSTRACT 

CNj ' Aims. We determine masses and mass-to-light ratios of five ultra-compact dwarf galaxies and one dwarf elliptical nucleus in the Fornax cluster 
,— H from high resolution spectroscopy. We examine whether they are consistent with pure stellar populations or whether dark matter is needed to 
explain their masses. 

' Methods. Velocity dispersions were derived from selected wavelength regions using a direct-fitting method. To estimate the masses of the 
^ C* | UCDs a new modelling program has been developed that allows a choice of different representations of the surface brightness profile (i.e. 
O j 1 Nuker, Sersic or King laws) and corrects the observed velocity dispersions for observational parameters (i.e. seeing, slit size). The derived 
dynamical masses are compared to those expected from stellar population models. 

Results. The observed velocity dispersions range between 22 and 30 kms~'. The resulting masses are between 1.8 and 9.5 x 10 7 M G . These, as 
well as the central and global projected velocity dispersions, were derived from the generalized King model which turned out to give the most 
^ stable results. The masses of two UCDs, that are best fitted by a two-component profile, were derived from a combined King+Sersic model. 
• • ■ The mass-to-light ratios of the Fornax UCDs range between 3 and 5 (M/Ly) e . The M/L v ratio of the dwarf elliptical nucleus is 2.5. These 
values are compatible with predictions from stellar population models. Within 1-2 half-mass radii dark matter is not dominating UCDs and the 
nucleus. An increasing dark matter contribution towards larger radii can not be ruled out with the present data. The M/L v ratios of some UCDs 
' » 1 ■ suggest they have intermediate age stellar populations. 
C\5 Conclusions. We show that the mass-to-light ratios of UCDs in Fornax are consistent with those expected for pure stellar populations. Thus 
UCDs seem to be the result of cluster formation processes within galaxies rather than being compact dark matter dominated substructures 
themselves. Whether UCDs gained their mass in super-star cluster complexes of mergers or in nuclear star cluster formation processes remains 
an open question. It appears, however, clear that star clusters more massive than about 5 x 10 6 M Q exhibit a more complex formation history 
than the less massive 'ordinary' globular clusters. 
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1. Introduction 

In the late 1990s a new type of very bright compact stellar sys- 
tem was discovered in the core of the nearby Fornax galaxy 
cluster (Hilker et al l 19991 Drinkwater et al. 120001 ). On ground- 
based images with typical resolution (~ 0.8 arcsec seeing) 
these objects look like stars since at Fornax distance globular 
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clusters (GCs) are unresolved. However, with absolute mag- 
nitudes in the range My = -13.5 to -12.3 they are about 2 
mag brighter than a> Centauri, the most massive globular clus- 
ter in the Milky Way, but about 3 mag fainter than the com- 
pact dwarf elliptical M32. Their luminosities are comparable 
to those of nuclei of dwarf ellipticals or late-type spirals (i.e. 
Lotz et al. 120041 Walcher et al. l2005l Cote et al. [20061 . Their 
sizes range between 10 and 30 pc in half-light radii (Drinkwater 
et al. |2003) which is about 10 times larger than ordinary GCs. 

In the last years, extensive surveys in the Fornax and Virgo 
cluster have revealed dozens of these compact objects (UCDs 
and bright GCs) in the magnitude range -13.5 < M v < -1 1.0 
(i.e. Mieske et al. l2002| and l2004l Hasegan et al. 120531 Jones et 
al. |2~006). They were dubbed "ultra-compact dwarf galaxies" 
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Table 1. Observation log and properties of observed objects. Vq is from Evstigneeva et al. (2006) and (V - I)q from Mieske et 
al. ( 1200 6) (except (V - 7)o for FCC303 which is from Karick et al. |20031 >. Both, Vo and (V - 7)o are corrected for foreground dust 
extinction (Schlegel et al. 119981 . The [Fe/H] values are from Mieske et al. (120061 1 . 



Name 


FCSS" name 


a(2000) 


(5(2000) 




(V-I) Q 


[Fe/H] 


Date 


Exp. time 


Seeing 


S/N* 


S/N* 






[h:m:s] 


[":':"] 


[mag] 


[mag] 


[dex] 




[s] 


["] 


4900A 5900A 


FCC303' 




3:45:14.10 


-36:56:12.0 


15.90 


1.09 




2000 Dec 20 


3 x 1004 


0.80 


3.6 


4.2 


UCD2 


J033806.3-352858 


3:38:06.33 


-35:28:58.8 


19.12 


1.14 


-0.90 


2000 Dec 23 


5 x 1800 


0.55 


4.5 


5.6 


UCD3 


J033854.1-353333 


3:38:54.10 


-35:33:33.6 


17.82 


1.22 


-0.52 


2000 Dec 20 


3 x 1000 


0.80 


4.7 


5.9 


UCD4 


J033935.9-352824 


3:39:35.95 


-35:28:24.5 


18.94 


1.13 


-0.85 


2000 Dec 24 


5 x 1800 


0.60 


4.8 


5.6 


UCD5 


J033952.5-350424 


3:39:52.58 


-35:04:24.1 


19.40 


1.00 




2001 Feb 12 


2 x 3285 


0.60 


4.4 


5.2 


UCD5 














2001 Feb 13 


2 x 3285 


0.55 


4.8 


5.8 


UCD5 














2001 Aug 19 


1 x 3285 


0.55 


4.6 


5.1 


UCD5 














2001 Sep 20 


1 X3119 


0.80 


5.3 


6.3 


UCD5 














2001 Nov 19 


1 x 3285 


0.85 


5.2 


6.4 



"Fornax Cluster Spectroscopic Survey (Drinkwater et al. 19991; ^measured for a single spectrum per resolution element (= 0.03); c Fornax 
Cluster Catalog (Ferguson [l989l > 



(UCDs) by Phillipps et al. (l200Tb . and "dwarf-globular tran- 
sition objects" (DGTOs) by Hasegan et al. (2005 ), expressing 
their uncertain origin. Various formation scenarios have been 
suggested to explain the origin of UCDs. The most promising 
are: 

1) UCDs are the remnant nuclei of galaxies that have been 
significantly stripped in the cluster environment (i.e. Bassino et 
al. [19941 Bekki et al. 120031 . This would mean that UCDs once 
resided in the centres of dark matter (DM) dominated cosmo- 
logical substructures. The isolated nuclei themselves, however, 
are probably not DM dominated since DM halos of galaxies 
can be stripped (Bekki et al. 2003 ). Moreover, present-day nu- 
clei of dwarf ellipticals are dominated by their stellar popula- 
tions, not by dark matter (Geha et al. 120021 ). 

2) UCDs formed from the agglomeration of many young, 
massive star clusters that were formed during merger events 
(i.e. Kroupa [19981 Fellhauer & Kroupa [20021). Hence, UCDs 
would be dark matter free stellar objects that were not created 
in the centres of galaxies. 

3) UCDs are the brightest globular clusters formed in the 
same GC formation event as their less massive counterparts 
(i.e. Mieske et al. 120021 Martini & Ho 2004]). Also in this case, 
UCDs would be objects possessing pure stellar populations. 

4) UCDs are the surviving counterparts of genuine compact 
dwarf galaxies that formed in small dark matter halos at the 
low mass end of cosmological sub-structure (i.e. Drinkwater et 
al. 2004). Nowadays UCDs survived the galaxy cluster forma- 
tion and evolution processes, e.g. they were not disrupted or 
accreted by larger galaxies until the present time. This would 
mean that they might still be dominated by dark matter. 

To distinguish the different formation scenarios one needs 
to study the detailed properties of UCDs, like structural pa- 
rameters, internal kinematics, together with ages and chemical 
compositions of their stellar population(s). In particular, the 
internal velocity dispersion cr of the stars within the UCDs 



is of interest, from which masses and mass-to-light (M/L) ra- 
tios can be estimated. A high virial M/L ratio would imply a 
compact cosmological progenitor rather than a merged stellar 
super-cluster. 

Recent high spatial resolution imaging and mid- to high- 
resolution spectroscopy shed some light on the properties of 
UCDs. First results on the velocity dispersion of UCDs in 
Fornax were presented in Drinkwater et al. (2003 ). These were 
based on a quick analysis of the same high resolution spec- 
tra that are the basis for this paper. The observed, projected 
velocity dispersions ranged between 24 and 37 kms -1 . Using 
the sizes of the UCDs - derived from high resolution sur- 
face brightness profiles - their masses and mass-to-light ratios 
(M/Ly) were roughly estimated. M/Ly values were of the order 
2-4 in solar units. This is slightly higher than the M/L ratio of 
globular clusters, but much lower than that of dwarf spheroidal 
galaxies of similar mass. In the present paper, we revise the 
formerly obtained masses and M/L ratios of the Fornax UCDs. 

In the Virgo cluster, Hasegan et al. ([2005 ) derived masses 
and M/L ratios for six DGTOs/UCDs. Three of those have 
relatively high M/Ly values between 6 and 9 at velocity dis- 
persions between 20 and 30 kms _I . This makes them the 
first UCDs with some evidence for dark matter, given their 
rather low metallicities (-1.5 <[Fe/H]< -0.7 dex). Fellhauer 
& Kroupa (120061) . however, point out that the tidal "heating" 
due to a close passage of an UCD through the central region of 
its host galaxy (pericenter between 100 and lOOOpc) may lead 
to an overestimation of its virial mass at all orbital radii due 
to the presence of unbound stars in the line of sight. The effect 
would occur mostly for UCDs with nearly radial orbits. Such 
orbits are indeed predicted for the remnant nuclei of stripped 
galaxies (Bekki et al. 120031 ) because of the selective way tidal 
stripping works. 

The structural parameters of the brightest Virgo and Fornax 
UCDs were presented in several recent papers (i.e. Hasegan 
et al. 120051 De Propris et al. 120051 Evstigneeva et al. l20"56l A 
striking result is that UCDs brighter than My = -11.5 follow 
a luminosity-size relation (Hasegan et al. 120051 Richtler[2006 
Kissler-Patig et al. 2006 Mieske et al. 2006 ), different from the 
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Table 2. Observation log and properties of standard stars that were used as velocity dispersion templates. They 
are ordered by their increasing metallicity. Most of the properties were taken from the SIMBAD database, see: 

|http : / / simbad . u-strsbg . f r/ Simbad] 



# 


Name 


Type" 


or(2000) 


(5(2000) 


V 


(B-V) 


[Fe/H] 




Date 


Exp. 


S/N* 


S/N 








[h:m:s] 


[":':"] 


[mag] 


[mag] 


[dex] 


[km s" 1 ] 




M 


4900A 


5900A 


1 


HD 21022 


G 


03:22:21.59 


-32:59:39.7 


9.20 


0.95 


-1.99 


110+3 


2000 Dec 20 


50 


50.0 


55.0 


2 


HD9051 


G 


01:28:46.47 


-24:20:25.3 


8.93 


0.82 


-1.50 


-73±1 


2000 Dec 20 


40 


52.0 


57.0 


3 


HD 41667 


G 


06:05:03.64 


-32:59:39.3 


8.53 


1.01 


-1.18 


302±10 


2000 Dec 23 


30 


49.0 


53.0 


4 


HD 20038 


G 


03:10:26.81 


-58:49:40.5 


8.90 


0.84 


-0.87 


27±10 


2000 Dec 24 


40 


52.0 


55.0 


5 


HD 17233 


SG 


02:43:52.03 


-54:47:28.3 


9.03 


0.79 


-0.60 


-18±10 


2001 Feb 13 


45 


49.0 


53.0 


6 


HR 296-1 


KOIII 


01:03:02.50 


-04:50:12.0 


5.43 


1.11 


-0.14 


15±? 


2001 Aug 19 


5 


58.0 


61.0 


7 


HR 296-2 
















2001 Sep 20 


5 


59.0 


62.0 


8 


HR 296-3 
















2001 Nov 19 


2.5 


53.0 


57.0 


"G: 


giant; SG: supergiant; ^measured for a single spectrum per resolution element (= 


0.03A) 
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Fig. 1. The six template spectra (see Table [2]i are shown in 
100A wide regions around the H/3 line (~ 4860A, left panels) 
and the Nal doublet (~ 5890A, right panels). The spectra are 
ordered by increasing metallicity as indicated. 
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Fig. 2. The five object spectra (see Table[TJi are shown in 100A 
wide regions around the H/3 line (~ 4885A, left panels) and the 
Nal doublet (~ 5920A, right panels). The 1-sigma error are 
plotted as well (grey lines). 



luminosity independent sizes of globular clusters (i.e. Jordan et 
al. 120051 ). Moreover, the light profiles of the brightest UCDs 
cannot be described by a simple King profile (Evstigneeva et 
al. [2006 ). Most of them show a cuspy profile in their centres 
and no tidal truncation within the observed surface brightness 
limit. Some of them exhibit small envelopes with effective radii 
upto 100 pc (e.g. Drinkwater et al. 120031 Richtler et al. l2005b . 

Concerning the chemical abundances of UCDs, recent 
work suggests that most of the bright Virgo UCDs are old 
and have metallicities ranging from [Z/H]= -1.35 to +0.35 
dex with an [a/Fe] value of ^0.3 (Evstigneeva et al. 2006), 
whereas the brightest Fornax UCDs are predominantly metal- 
rich ([Fe/H]^ -0.5 + 0.5 dex; Mieske et al. 12006] ). Their ages 
and a-abundances are not constrained very well. Future surveys 



will clarify whether some UCDs may exhibit intermediate-age 
stellar populations. 

In this paper, we determine dynamical masses of five of 
the brightest UCDs and one nucleus of a dwarf elliptical in the 
Fornax cluster. The work is based on the detailed analysis of 
high resolution spectroscopy that was obtained with the Very 
Large Telescope Ultraviolet and Visual Echelle Spectrograph 
(VLT/UVES) at Paranal/Chile (and Keck spectroscopy for one 
of the UCDs). The derivation of structural parameters for the 
UCDs - which are necessary to estimate virial masses - are 
taken from Evstigneeva et al. (2006 ). In section 2 of the present 
paper, the observations and data processing are described. 
Section 3 is dedicated to the velocity dispersion measurements, 
and section 4 to the derivation of masses and mass-to-light 
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Fig. 3. The "blue" part (upper panel) and the "red" part (lower panel) of the spectrum of UCD4 are shown. Dark regions mark 
the eight direct fitting regions for velocity dispersion measurements (see Table |3}. The grey wavy lines are the 1-sigma error 
spectra. 



ratios by dynamical modelling and their comparison with ex- 
pected values from stellar population models. Finally, in sec- 
tion 5 the main results are summarized. Throughout the paper 
we use a distance modulus of 31.39 mag for the Fornax cluster 
(Freedmanetal. l200TT) . 

2. Observations and data reduction 

2.1. Observations 

The observations were performed between December 2000 and 
November 2001 with the VLT/UT2 at Paranal (ESO, Chile) in 
service mode. The instrument in use was the Ultraviolet and 
Visual Echelle Spectrograph (UVES), and an array of two at- 
tached 2x4 k CCD chips. The data were read out with a 2x2 
binning, resulting in a spatial scale of 0.36"/pixel. The setting 
used (RED arm, /t ce ntrai = 5200A) provided a wavelength cov- 
erage of 4200-6200A, with a gap between 5160 and 5 240 A 
(blind region beween the two CCD chips). The slit was chosen 
to have a width of 1 "which gives a resolving power of 37000- 
45000. The final resolution ranged between 2.4 and 3.1 pixels 
(or 0. 12 and 0. 16A). It depends on the seeing which varied be- 
tween 075 and 079. 

In total, four ultra-compact dwarfs in the magnitude range 
17.8 < V < 19.4 mag and one nucleus of a dwarf ellipti- 
cal (V nuc - 19.5) were observed. Most of the observations 
were taken in three nights in December 2000. Only for one ob- 
ject (UCD5), the single exposures were taken in five different 
nights between February and November 2001. The total inte- 
gration time for the UCDs varied between 50 minutes and 6 
hours. Table Q] gives a short logbook of the observations, to- 
gether with the achieved S/N at two different wavelengths (per 
0.05A pixel). Also the properties of the observed objects are 
summarized in this table. 

In the same run, spectra of six standard stars were taken 
with the same instrument configuration. They mostly are gi- 
ants and cover a metallicity range of -2.0 <[Fe/H]< -0.1 dex. 



Table 3. Wavelength regions for direct fitting of velocity dis- 
persions and prominent lines therein. The values are given in 
the restframe. 



Id 


^start 


-ieiid 


Absorption lines 




[A] 


[A] 




1 


4523.2 


4632.7 


Mgl 45 71 


2 


4816.8 


4931.3 


CIV, H/? 


3 


4946.2 


5075.6 


Fe 504 i 


4 


5274.6 


5384.1 


Fe 53 7i 


5 


5583.1 


5682.6 




6 


5866.8 


5916.5 


Nal (D2), Nal (Dl) 


7 


5921.5 


6031.0 




8 


6090.7 


6155.4 


Call, CaI2 



Their properties are summarized in Table [2] and excerpts of 
their spectra around the Hfi line and the Na doublet are shown 
in Fig. Q] The same spectral regions (however redshifted) for 
the UCDs and FCC303 are shown in Fig. El 

2.2. Data reduction 

The European Southern Observatory (ESO) provides pipeline 
reduced data for UVES echelle spectroscopy. The UVES con- 
text in the MIDAS environment is used for that reduction. Due 
the low S/N of our data, the pipeline reduced data with stan- 
dard settings can not be trusted right away. Therefore, we de- 
cided to re-reduce these data with the same pipeline, carefully 
controlling all steps and optimizing the setup parameters. Also, 
we implemented a different more reliable cosmic ray removal 
routine as the one used by the pipeline. 

The first steps, namely first guess solution, defining of order 
positions, wavelength calibration, were performed with the de- 
fault settings of the UVES pipeline. Then, before the order trac- 
ing and extraction of the spectra, the cosmic rays in the low S/N 
science exposures were rejected by the LACOSMIC (Laplacian 
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Fig. 4. The direct fitting regions for the spectrum of UCD4 are 
shown (grey lines) together with the best fitting template (black 
lines). The name of the template, the velocity dispersion of the 
gaussian kernel and the number of the direct fitting region are 
indicated in each panel. 



Cosmic Ray Identification) routine from van Dokkum (2001) 
which worked very well. This was necessary to avoid jumps 
in the tracing as a result of cosmic rays close to the spectra. 
For the further reduction, the inter-order background subtrac- 
tion method MEDIAN, the OPTIMAL extraction, flatfield cor- 
rection in the extracted pixel-order space, sky subtration, and 
merging of overlapping orders were chosen. Also an error spec- 
trum containing the standard deviation of the reduced spectra 
was extracted for each object. 

Since some spectra of the same object were observed at 
very different dates, a correction for the heliocentric velocity 
was necessary before combination. For UCD5, the heliocen- 
tric corrections differed by 35.08 km s _1 between individual 
exposures. Since the instrumental spectral resolution is about 
0.03 A/pixel, a noticable shift in wavelength (10% = 0.003 A) 
already occurs for a velocity difference of > 0.15 km s . 
The combination was done with the SCOMBINE under the 
NOAO.IMRED.ECHELLE package of IRAF. Prior to averaging 
the single spectra they were scaled by their mode, and pixel 
values deviating more than 3 sigma from the average were re- 
jected. In the same step the spectra were rebinned to a resolu- 
tion of 0.05 A/pixel. The averaging and rebinning increased the 
S/N of a resolution element to 11-19 for the UCDs, depending 
on the object and the wavelength (see Table |4). 

The standard star spectra have a high S/N (~ 50-60) and 
were reduced with the standard settings of the pipeline. They 
were rebinned to the same resolution as the science spectra. 
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Fig. 5. The Na doublet region for the five studied objects is 
shown (UCD2, UCD3, UCD4, UCD5 and FCC303 from top 
to bottom). The bold lines are averages of the three best-fitting 
templates for this region (region 6, see Table [3]). The corre- 
sponding velocity dispersions are given on the left. Thin lines 
represent the averages of the three best-fitting templates in the 
regions 4, 5 and 7. Their velocity dispersions are given on the 
right. 

3. Velocity dispersion measurements 

The kinematic analysis of the spectra was performed using 
van der Marel's direct-fitting method (van der Marel & Franx 
1993jl. First, the spectra were placed on a logarithmic wave- 
length scale and normalized. Then, the standard star spectra 
were convolved with Gaussian velocity dispersion profiles in 
the range 2 to 60 km s and a step size of 1 km s _1 . 

For the comparison of the templates with the UCD spectra, 
wavelength regions with prominent absorption features were 
selected. The restframe wavelengths of the direct fitting regions 
and some prominent lines they cover are summarized in Table[3] 
and are shown in Fig. [3] Additionally, direct fitting on a combi- 
nation of three selected regions was performed for the bluer 
part (regions 1, 2 and 3) as well as for the redder part (re- 
gions, 4, 5 and 7) of the spectrum. All object spectra (UCDs 
and FCC 303) were fitted with all sets of smoothed template 
spectra (see Table|2]i. The best-fitting Gaussian velocity disper- 
sion profile is determined by x 2 minimization in pixel space. 
As an example, the best fits to the spectrum of UCD4 for the 
eight direct fitting regions are shown in Fig. [4] In Table [4] the 
results for the three best-fitting templates (average values) are 
given for all fitting regions. The identification numbers of the 
three templates are given in brackets. In most cases, the de- 
rived velocity dispersions of these templates are consistent with 
each other within +2 km s _1 . The template star HR296 was 
observed in three different nights. The velocity dispersions as 
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Table 4. Results of velocity dispersion measurements in all direct fitting regions (same identification numbers as in Table|3j. On 
the left, the results from the blue part of the spectrum are shown, and on the right the ones from the red part. The identification 
numbers of the three best fitting templates are given in brackets below each cr v value. The S/N are those for the combined spectra 
at the given wavelength. 



Name 


S/N 




cr v 


[km s- 1 ] 




S/N 






(T v [km s 1 ] 








4900A 


1 


2 


3 


1+2+3 


5900A 


4 


5 


6 


7 


8 


4+5+7 


UCD2 


14.3 


24.4±1.5 


24.8+1.5 


24.7+1.4 


25.3+0.8 


15.9 


22.5+1.7 


22.8+2.3 


26.5+1.7 


24.3+3.0 


19.4+2.0 


24.4+1.2 






(3,4,5) 


(5,4,3) 


(8,4,3) 


(3,4,5) 




(6,5,4) 


(3,5,8) 


(8,2,5) 


(6,4,3) 


(8,4,5) 


(8,4,3) 


UCD3 


11.3 


23.4+1.2 


25.8+1.1 


26.1 + 1.0 


25.4+0.7 


13.0 


29.2+1.4 


29.6+2.0 


44.3+1.6 


26.8+3.0 


32.3+2.6 


30.0+1.1 






(8,3,4) 


(5,7,4) 


(8,5,4) 


(8,5,4) 




(7,5,3) 


(5,8,3) 


(2,8,5) 


(6,3,4) 


(8,5,4) 


(7,3,4) 


UCD4 


14.9 


27.6±1.9 


30.1 + 1.7 


27.5+1.6 


29.3+1.0 


15.9 


29.4+2.1 


26.4+2.1 


30.2+1.8 


29.2+7.9 


31.1+3.2 


31.3+1.6 






(8,4,3) 


(4,5,3) 


(4,5,2) 


(3,4,5) 




(3,5,8) 


(3,2,1) 


(2,5,8) 


(6,5) 


(5,3,4) 


(3,8,4) 


UCD5 


17.6 


24.1 + 1.8 


26.7+1.9 


23.5+1.5 


25.4+1.0 


19.1 


20.1 + 1.7 


19.0+1.8 


23.8+1.9 


20.3+2.7 


23.1+4.6 


20.7+1.1 






(3,4,5) 


(3,4,2) 


(3,4,2) 


(3,4,5) 




(3,4,1) 


(4,5,3) 


(6,1,2) 


(7,4,3) 


(4,8,5) 


(4,3,8) 


FCC303 


7.1 


19.9+5.3 


24.7+2.8 


26.8+2.2 


29.6+1.8 


9.3 


24.4+2.6 


23.3+2.7 


34.0+4.0 


6.5+3.9 


31.6+6.8 


28.1+2.1 






(3,2,1) 


(4,2,5) 


(2,4,5) 


(3,4,2) 




(4,8,5) 


(3,5,7) 


(4,5,8) 


(3,2,6) 


(5,4,3) 


(3,4,2) 



derived from the three individual velocity dispersion fits with 
this star agree within +1 km s -1 with each other in all cases. 
Thus no systematic effects larger than ±2 km s -1 either due to 
the use of different templates or observations in different nights 
are expected. The internal velocity dispersions of the stars are 
expected to be below 3 km s~' , and thus do not affect the veloc- 
ity dispersion measurements. Note that one resolution element 
of 0.05 A/pixel corresponds to 2.5-3.5 km s -1 . 

The best-fitting templates are in most cases the ones with 
metallicities between -1 and dex. This should be expected 
since line index measurements of the Fornax UCDs 2, 3 and 4 
revealed iron abundances of [Fe/H]= -0.90, -0.52 and -0.85 
dex, respectively (Mieske et al. (2006). 

The scatter of velocity dispersion measurements between 
the different direct fitting regions is very different for the bluer 
and redder part of the spectrum. Whereas in the bluer part the 
scatter is of the order +1.5 kms -1 , the regions 6 and 8 of the 
redder part, which include the Na doublet and a Ca line, reveal 
large discrepancies towards high cr v values. This effect is illus- 
trated in Fig. [5] where the best fit to the Na doublet region is 
compared to the expected fit in this region when taking the av- 
erage <j v of regions 4, 5 and 7. This discrepancy cannot easily 
be explained. Probably abundance mismatches and/or stellar 
population mixes play a role. The c-abundances of UCDs in 
Virgo are known to have super-solar values (Evstigneeva et al. 
2006). One might also speculate whether the extra absorption 
might be caused by intra-cluster material in the line of sight. 
Future studies shall clarify this issue. 

Table|4]also shows that there seems to exist a slight system- 
atic effect between the velocity dispersion measurements in the 
bluer and redder part of the spectrum (excluding regions 6 and 
8) for some objects. UCD3 has systematically higher cr values 
in the red part as compared to the blue part, whereas for UCD5 
the contrary is the case. Also here, template mismatches and 
contribution of stars from different evolutionary stages might 
be the reason. 

To define a final velocity dispersion for each object we con- 
sider three cases of averaging the individual values: 1) simply 
the mean of all fits; 2) the mean of all fits but with a 2-sigma 



Table 5. Average values for radial velocity and velocity dis- 
persion measurements: (1) - mean of all fits; (2) - as (1), but 
with 2-sigma clipping; (3) - mean excluding Na and Ca region 
and combined region. The values in bold face are those used for 
further analyses. The velocity dispersion of UCD1 was derived 
by Drinkwater et al. (2003) from Keck spectroscopy. 



Name Vheiio o"v °"v o~ y 

[kms -1 ] [kms -1 ] [kms -1 ] [kms -1 ] 
(1) (2) (3) 



UCD2 


1205.7+0.8 


24.4+1.6 


24.7+1.0 


24.1 + 1.0 


UCD3 


1448.3+1.5 


27.6+4.9 


26.6+2.3 


26.2+2.2 


UCD4 


1839.8+3.8 


29.3+1.7 


29.2+1.5 


28.4+2.0 


UCD5 


1245.4+2.1 


22.9+2.6 


22.9+2.6 


22.4+2.9 


FCC303 


1928.6+3.1 


26.1+5.6 


27.1+3.1 


23.1+6.0 


UCD1 






32.0+1.0 





clipping of outliers; 3) mean excluding the Na and Ca region as 
well as the combined regions. The results are listed in Table [5] 
They are consistent with each other within the errors. We adopt 
the average velocity dispersion of case 2) as final value for fur- 
ther analyses. 

For UCD1 no UVES spectroscopy is available. Yet its ob- 
served velocity dispersion is known from Keck spectroscopy 
(Drinkwater et al. 2003 ). It was derived from the CaT region 
using a cross-correlation method (Tonry & Davis 119791 as 
implemented in RVSAO/IRAF) and only one stellar template 
(G6/G8IIw type). The determined value is given in Table [5] 
and is used for the mass modelling (see next section). Note 
that in Drinkwater et al. (2003) also the velocity dispersions of 
UCD2 to UCD5 and FCC303 were determined via the cross- 
correlation method. Since in this preliminary study only the Na 
doublet region was used the obtained velocity dispersions are 
2- 10 km s -1 higher than the adopted values of the present study 
(see discussion on the Na doublet above). 



M. Hilker et al.: Dynamical masses of of ultra-compact dwarf galaxies in Fornax 



7 



4. Mass estimates and mass-to-light ratios 

The masses and mass-to-light ratios of the UCDs are important 
physical parameters for the understanding of their origin. In 
particular, the mass-to-light ratio (M/L) is an indicator for pos- 
sibly existing dark matter and/or violation of dynamical equi- 
librium. If UCDs were the counterparts of globular clusters 
- thus a single stellar population without significant amounts 
of dark matter - one would expect M/L values as predicted 
by standard stellar population models (e.g. Bruzual & Chariot 
2003 Maraston 2005 ). If UCDs were of cosmological origin - 
thus formed in small, compact dark matter halos - they might 
still be dominated by dark matter and show a high M/L value. 
Mass-to-light ratios that are larger than expected from single 
stellar populations can, however, also be caused by objects that 
are out of dynamical equilibrium, e.g. tidally disturbed stellar 
systems (Fellhauer & Rroupa l2006t or by stellar populations 
that formed with an unusual initial mass function. 

In the following section we describe how the masses and 
mass-to-light ratios were determined. 

4. 1. Mass modelling 

The ingredients for the mass modelling of UCDs are their mea- 
sured velocity dispersions and their observed structural param- 
eters. The light profiles from Hubble Space Telescope imag- 
ing of several UCDs in Virgo and Fornax are presented in 
Evstigneeva et al. (2006). It has been shown that a simple King 
profile often is not the best choice to represent their surface 
brightness profiles. Instead, Nuker laws or generalized King 
profiles (with variable a) fit the light profiles much better (see 
formulae in Evstigneeva et al.). Most of the mass estimators 
available in the literature, however, are based on the assump- 
tion of a King model with exponent a — 2. 

Not to be restricted only to King profiles a general approach 
has been developed. It calculates an objects' mass considering 
various representations of their light profiles. Thus, the follow- 
ing steps are performed: 

1 . Fitting the observed luminosity profile by a given density law 
(Nuker, King, generalized King, Sersic or a two-component 
King+Sersic function). 

2. Deprojection of the fitted 2-dimensional surface density pro- 
file by means of Abel's integral equation into a 3-dimensional 
density profile. 

3. Calculation of the cumulated mass function M(< r) and the 
potential energy (p(r) from the 3-dimensional density profile. 
From these the energy distribution function f(E) is then cal- 
culated with the help of eq. 4- 140a from Binney & Tremaine 
( 1987 ), assuming isotropic orbits for the stars. 

4. Creation of an iV-body representation of the UCD using the 
de-projected density profile and the distribution function. For 
every model, 100.000 test particles were distributed. 

The modelling is based on the assumptions of spherical 
symmetry, an underlying isotropic velocity distribution and 
no bias between mass and light (i.e. mass segregation is ne- 
glected). 

Spherical symmetry is justified since most UCDs are close 
to round with projected ellipticities in the range 0.01 < e < 



0.24 (see Table 6, e — 1 - b/a). Their ellipticity distribution 
is consistant with that of Milky Way and NGC5128 globular 
clusters (Evstigneeva et al. 2006). Concerning the modelling, 
the potential is always more spherical than the surface density, 
so a small ellipticity does not have much influence on the mea- 
sured velocity dispersion. 

The assumption of isotropy is justified by recent proper mo- 
tion and radial velocity observations of a> Centauri, which show 
that the velocity dispersion of this cluster is close to isotropic 
inside two half-mass radii (van de Ven et al. 120061 . Since oj Cen 
might be regarded as a Local Group analogue for a UCD, the 
UCDs in our sample could have isotropic velocity dispersions 
as well. In addition, two-body relaxation will have isotropised 
the orbits of stars, at least in the inner parts of the studied 
UCDs. 

The assumption of mass following light is a critical issue. 
Whereas mass segregation of the stars can be neglected be- 
cause of the long relaxation time of these massive objects, an 
extended dark matter (DM) halo can not be excluded a priori. 
From the modelling point of view the addition of a DM halo is 
straightforward. Having, however, only one measured velocity 
dispersion (instead of a velocity dispersion profile), there is a 
large freedom in the scale length and total mass of the DM halo 
we model to match the observed velocity dispersion as well as 
the stellar density profile. In the cluster outskirts one could add 
nearly any amount of DM without changing the central veloc- 
ity dispersion significantly. Not much can be learned from such 
an exercise. Therefore, we postpone the modelling of UCDs 
with extended dark matter halos to the future when spatially 
resolved velocity dispersion profiles become available. 

Besides the projected profile parameters, the by us devel- 
oped program expects as input the total mass of the stellar sys- 
tem and the number of test particles, and creates a list of x, 
y and z positions and v x , v y and v, velocities for all particles 
that correspond to the specified structural parameters and the 
given mass. From this output file the central as well as global 
projected velocity dispersion can be calculated. The modelled 
projected half-mass radius, r^p is given as well for comparison 
with the observed value. In Fig. [6] the 3-dimensional density 
distribution, the cumulative mass distribution and the projected 
velocity dispersion profile is shown in for four different model 
representations of UCD4 (see Table |5J. The same is shown for 
the Nuker and King+Sersic models for UCD3 in Fig. [7] 

In case of the Nuker and Sersic functions, the models were 
truncated at large radii to avoid the unphysical infinite exten- 
sions of UCD light profiles. The truncation radius of the Nuker 
model was fixed at 2 kpc. Depending on the steepness of the 
outer power law the total mass estimate can depend critically 
on the exact value of the truncation radius. See, for example the 
rising cumulative mass of UCD4 in Fig. [6] (middle panel). The 
truncation radius of the Sersic model was set to 20 times the 
effective radius, thus ranging between some hundred parsecs 
and a few kiloparsecs for the UCDs in our sample. In some 
cases, the adopted truncation radius might be too short leading 
to lower values of the total mass as compared to other models 
(see Fig. [6) middle panel). 

The true tidal radii of the UCDs depend on their distances 
to the cluster center Rq and the mass m c of the cluster (UCD). 
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12 3 



log(radius) [pc] 

Fig. 6. The model output of four different light profile repre- 
sentations of UCD4 is shown; the line types are: black solid 

- Nuker law, black dot-dashed - Sersic law, grey long dashed 

- King profile, grey short dashed - generalized King profile. 
The upper panel shows the 3-dimensional density distribution, 
the middle panel the cumulative mass distribution, and the 
lower panel the projected velocity dispersion profile of the z- 
component of the 3-dimensional velocity. The vertical short 
dashed lines indicate the radii of 0.5 (left) and 1 pixel. The 
dotted line is half the slit width, and the bold dashed line shows 
the projected half-light radius of UCD4. 



UCD3 




12 3 

log(radius) [pc] 

Fig. 7. The model output of two different light profile repre- 
sentations of UCD3 is shown; the line types are: solid - 2- 
component King+Sersic profile, dashed - Nuker law. The upper 
panel shows the 3-dimensional density distribution, the middle 
panel the cumulative mass distribution, and the lower panel the 
projected velocity dispersion profile of the z-component of the 
3-dimensional velocity. The vertical short dashed lines indicate 
the radii of 0.5 (left) and 1 pixel. The dotted line is half the slit 
width, and the bold dashed line shows the projected half-light 
radius of UCD3. 



It can be estimated by the formula: r, = (Gm c /2v 2 . )^ 3 /??/ 3 , 
where Vdrc is the circular velocity of the cluster potential and G 
the gravitational constant. For the Fornax cluster, the circular 
velocity is about 420 kms -1 (e.g. Richtler et al. 2004), and the 
projected distances of the UCDs are between 30 and 180 kpc. 
The resulting lower limits (because of projected distances) of 
the tidal radii range between 1 and 2 kpc, agreeing with the 
adopted truncation radii. 

After the creation of the UCD model the velocity disper- 
sion as seen by an observer is calculated. To achieve this, the 
following steps are performed: 



1 . All test particles are convolved with a Gaussian whose full- 
width at half-maximum (FWHM) corresponds to the observed 
seeing. 

2. The fraction of the 'light' (Gaussian) that falls into the slit 
at the projected distance of the observed object (the size of the 
slit in arcsec and the distance to the object in Mpc are input 
parameters) is calculated. 

3. These fractions are used as weighting factors for the veloci- 
ties. All weighted velocities that fall into the slit region are then 
used to calculate the 'mimicked' observed velocity dispersion 

Of course, the total 'true' mass of the modelled object, 
M tlue , that corresponds to the observed velocity dispersion, cr t, s 
is not known a priori. One has to start with a first guess of the 




I I I i r 

UCD2 (gen. King) 



U 

Jb c - 



-HOC 



-40G 
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Fig. 8. The x-y-plane of the modelled generalized King rep- 
resentation of UCD2 is shown. Grey dots are all 100.000 test 
particles. Black dots are those 'stars' whose centres fall into 
the analysed slit area. The grey circle indicates the projected 
half-light radius of UCD2. The Gaussian in the lower right cor- 
ner corresponds to the convolution kernel used to estimate the 
fraction of light that falls into the slit region for each 'star'. 



-100 -50 50 100 

v z, model [km/s] 

Fig. 9. Velocity histograms of the z-component velocities for 
the Nuker (top panel) and generalized King (lower panel) mod- 
els of UCD4. The bin width is 1 kms" 1 . Black solid lines in- 
dicate the velocity distributions of all test particles, the light 
grey lines those of the 'stars' with centres in the slit area, and 
the grey lines the luminosity weighted number counts of those 
'stars' whose 'light' falls into the slit area. The corresponding 
velocity dispersions are given in the plots. 



total mass, M guess , from which the 'true' mass can be calculated 

as Mime = M gU e S s • (o" bs/cr mod ) 2 . 

As an example, in Fig.[8]the x-y-plane of the modelled gen- 
eralized King representation of UCD2 is shown. The black dots 
indicate 'stars' whose centres fall into the analysed slit area 
(1x3 arcsec); the Gaussian represents the point spread function 
profile of the 'stars' caused by the seeing conditions during the 
observations (see TableQ}. In Fig.|9]the velocity histograms of 
all particles for two model representations of UCD4 are shown. 
The particles in the histograms are the ones in the slit area 
and the luminosity weighted number counts of those particles 
whose 'light' falls into the slit area. One can see that the seeing 
hardly changes the velocity dispersion for stars within the slit 
region. This is because most of the objects' light is contained in 
the analysis area (e.g. the half-light radii of UCDs are smaller 
than the slit width). Though, the total as well as central veloc- 
ity dispersion (see also Table [6]l can differ significantly from 
the observed one due to the extended light profiles and central 
concentration of the models, respectively. 

4.2. Dynamical masses 

The results of the modelled cluster masses and velocity 
dispersions for King, generalized King, Sersic, Nuker and 
King+Sersic functions are presented in Table [6] together with 
the surface brightness profile model parameters which were 
taken from Table 9 in Evstigneeva et al. (120061 1. 



The errors of the modelled masses were estimated from the 
uncertainty in the observed velocity dispersion, Act and from 
typical uncertainties of model parameters for surface bright- 
ness profiles (see Evstigneeva et al. [2006). On the one hand, 
maximum and minimum model masses were estimated us- 
ing the minimum and maximum observed velocity dispersion 
(c"obs+Acr) and (cr b s -Acr), respectively. The average of the dif- 
ferences (Afmax -Mtrue) and (M true -M m i n ) defines the first mass 
error. On the other hand, models were created that simulate 
o"obs from the profile parameters plus/minus their uncertainties. 
The maximum and minimum mass deviations from M tlue de- 
fine the second error. Both errors then were summed to derive 
the total mass error. 

The mass-to-light ratio errors were propagated from the 
mass errors and the uncertainty in the luminosity (assumed 
to be 0.05 mag in the absolute magnitude). The luminosities 
were derived from the apparent, extinction corrected V mag- 
nitude as given in Table [TJ a distance modulus to Fornax of 
(m - M)v — 31.39 mag (Freedman et al. 200J} and a solar 
absolute V magnitude of My >@ = 4.85 mag. 

The error of the central velocity dispersion was estimated 
from the observational error plus the scatter of modelled ve- 
locity dispersions in annuli of 0.5 parsec within the central 5 
pc for each object. The error of the global velocity dispersion 
is the sum of the observational error and the uncertainties as 
propagated from the mass error of the profile parameters. 
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Table 6. Mass modelling results for different light profile representations. The modelled projected half-light radius r^ p , central 
velocity dispersion cr , global velocity dispersion cr, mass and mass-to-light ratio in the V-band are given. The adopted surface 
brightness profile parameters are taken from Evstigneeva et al. (2006 ). Radii (R/,, R e , R c and R t ) are in pc; //o.King an d /^.Sers. are 
in mag arcsec~ 2 in the V band. In the last rows the results for the virial mass estimator are given. The effective radius R e g was 
taken from Table 10 of Evstigneeva et al. ( 1200 6 ), the global velocity dispersion from the generalized King models and King + 
Sersic models in two-component fits case. Also the ellipticies (e = 1 - b/a) from Evstigneeva et al. are shown for reference. The 
values in bold face were adopted for further analyses. 
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UCD1 


36.7 


0.34 


3.03 


0.91 






39.1 ±3.0 


51.4 ±3.0 


25.2 ± 2.0 


2.59 ±0.32 


4.03 ± 0.53 


UCD2 


4.3 


0.73 


2.38 


0.58 






39.4 ±3.0 


33.6 ± 1.0 


20.4 ± 2.0 


2.39 ±0.38 


3.45 ± 0.59 


UCD3 
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28.3 ± 1.5 
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2.40 ± 0.58 


3.47 ± 0.86 


UCD4 
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c 
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UCD1 


1.8 
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2.32 
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22.4 ± 1.0 
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UCD2 
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4.66 ± 1.32 
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0.05 




29.5 


25.0 ±3.4 
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0.24 




31.2 


18.7 ±3.2 
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0.03 




660.0 


16.5 ±3.9 
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The masses and M/L values of the different models in gen- 
eral agree with each other within the errors. In Fig. \W[ the 
luminosity versus mass plane is shown for all model results 
and the virial mass estimate (see below). The mass-to-light ra- 
tios of most UCDs scatter around M/Ly = 4. On average, the 
masses of the Sersic profile are slightly lower than those de- 
rived from the other profiles, mainly caused by the adopted 
truncation radius (see discussion above). The masses of the 
Nuker profile are sensible to the slope of the outer power law 
and the adopted truncation radius. The shallow outer slope of 
UCD4 causes a high mass estimate, whereas the rather steep 
outer slope of UCD1 results in a low mass estimate as com- 
pared to the other models. Besides the uncertainties in the trun- 
cation radius, Nuker and Sersic mass estimates suffer from the 
uncertainties in the central parts of the profile. The cuspy nature 
of the density profile depends on the light distribution within 



the central 2-4 pixels, an area for which it is difficult to make a 
model independent of the assumed fitting function (see Fig. 3 in 
Evstigneeva et al. 2006). Most sensitively, the central velocity 
dispersion ctq depends on the model shape within the innermost 
pixels. The derived ctq values for Nuker or Sersic models are 
3-10 km s higher than those for the King models. This has to 
be taken into account when using central velocity dispersions 
in fundamental plane plots. 

The most stable results for one-component profile fits are 
those of the generalized King model because they do not suffer 
from the uncertain profile shape in the inner pixels (like Nuker 
and Sersic) and fit the outer regions much better than the 'nor- 
mal' King models. Concerning two-component profiles, King 
core+Sersic halo models are the most reliable. Therefore, the fit 
values of these models were adopted for further analyses (see 
numbers in bold face in Table |6]). 
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jq L Bruzual & Chariot (2003)_ 



log(L v ) [L s ] 

Fig. 10. Luminosity vs. mass plot for five Fornax UCDs (the 
data for UCD1 were taken from Drinkwater et al. (2003). 
The symbols indicate the following models: grey open cir- 
cles - King, filled circles - generalized King, grey triangles - 
Sersic; grey squares - Nuker, open hexagons - King+Sersic, and 
crosses - virial mass estimator. Lines of constant mass-to-light 
ratio are shown as indicated. 



The results of our dynamical modelling were also com- 
pared with the virial mass estimator by Spitzer ( 1987 ), M v ; r « 
9.75(fl e ffcr 2 )/G, where cr is the global projected velocity dis- 
persion and R e g the half-light radius. This method also assumes 
a constant Mj L ratio as function of radius and an isotropic ve- 
locity distribution. The results are shown in Table [6] The errors 
in the virial mass estimates were propagated from the errors in 
/? e ff and cr. The virial and modelled masses agree within the 
errors. 



4.3. Comparison with masses from stellar 
population models 

In turn the masses and mass-to-light ratios as derived from dy- 
namical modelling are compared with the M/L values expected 
from stellar population models in order to see whether UCDs 
can be explained by pure stellar systems or whether they would 
need dark matter within the central 1-3 half-mass radii to ac- 
count for missing mass. 

There are stellar population models available from vari- 
ous groups. Here we present the comparison with Bruzual & 
Chariot (2003 ) and Maraston (120051 ) models. Mass-to-light ra- 
tios are given for single stellar populations (SSPs) of various 
ages and metallicities. The metallicity in the Maraston mod- 
els is a total metallicity [Z/H] including possible [a/Fe] over- 
abundances. [Z/H] and the iron abundance [Fe/H] are related as 




jq ^Maraston (2005) 
I [<x/Fe] = 0.3 




jq Z_ Maraston (2005) 
I [a/Fe] = 0.3 




_L 



I 1 1 1 I 1 1 
. Chabrier IMF 



Kroupa IMF 



Salpeter IMF 




_L 



-2 -1 
[Z/H] [dex] 



1 1.2 1.4 

(V-I) [mag] 



Fig. 11. The dynamical M/Ly ratios of Fornax UCDs (filled 
circles) and the nucleus of FCC303 (filled triangle) are com- 
pared to expected M/Ly values from stellar population mod- 
els. Also shown are data for Virgo UCDs (open squares from 
Evstigneeva et al. (2006). In the left panels M/Ly is plotted 
versus metallicity [Z/H], in the right panels versus the colour 
(V — T)o. The references to the models and the adopted IMF are 
indicated in the panels. The three lines in each plot represent 
ages of 15, 13 and 5 Gyr from top to botton, respectively. For 
the Fornax UCDs an a abundance of [or/Fe]= 0.3 was adopted 
to derive their [Z/H] values. Note, however, that their a abun- 
dance might be different (solar) if they possess a significant 
contribution of intermediate age stellar populations. 



following: [Fe/H] = [Z/H] -0.94-[a/Fe] (Thomas et al. 120031 . 
[Z/H] is in agreement with the Zinn & West ( 1984) scale. 

In Fig. QT| the dynamical masses of Fornax UCDs and the 
nucleus of FCC 303 are compared to stellar population mod- 
els in the [Z/H] versus M/Ly and (V - I) versus M/Ly plane. 
For comparison we also include data for seven of the brightest 
Virgo UCDs (Evstigneeva et al. 2006). Three sets of models 
with different assumption for the initial mass function (IMF) 
are presented. The Bruzual & Chariot models with a Chabrier 
IMF, and the Maraston models with a Kroupa as well as a 
Salpeter IMF. These IMFs were chosen to illustrate the wide 
range of predictions from various SSP models. The colours of 
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the Virgo UCDs were taken from Evstigneeva et al. (2006). The 
ages, metallicities and a abundances were estimated from diag- 
nostic plots of line indices in Evstigneeva et al. (Figs. 11 & 12). 
Most Virgo UCDs are old (8 to 15 Gyr), and most of them show 
super-solar [a/Fe] values around 0.3 dex. The iron abundances 
and colours of the Fornax UCDs were taken from Mieske et 
al. d20061 l. No age and or-abundance estimates are published 
for them. The relatively high Balmer line H/3 indices of UCD3 
and UCD4 (see Table 1 in Mieske et al. 120061 ) might point to 
a contribution of young to intermediate age stellar populations 
within these UCDs. Note that this statement has to be taken 
with caution since no proper calibration of the line indices was 
performed in Mieske at al.'s paper. 

Fig. [TTI shows that the dynamical masses of all UCDs and 
the nucleus can be reconciled with stellar population models if 
taking the 'right' set of models. The Bruzual & Chariot mod- 
els with a Chabrier IMF tend to predict MjLy values that are 
lower than the measured ones. Taking these models at face 
value one might conclude that most UCDs need some amount 
of dark matter to explain their higher MjLy values. When tak- 
ing Maraston models with a Salpeter IMF, however, even the 
extreme case of Strom417 (a low mass UCD or massive GC) 
with a dynamical MjLy value around 6 (Evstigneeva et al. 
2006 Hasegan et al. [2005 > can be explained by an old stellar 
population. If these models are correct for all UCDs, many of 
them are best fitted by intermediate age (5-8 Gyr) stellar popu- 
lations. In the Fornax cluster, UCD2, UCD3 and the nucleus of 
FCC303 would be candidates for young ages, whereas UCD1, 
UCD4 and UCD5 still are compatible with old stellar popula- 
tions. In the Virgo cluster, VUCD3, VUCD5 and VUCD6 are 
possibly younger than the other UCDs. This, however, contra- 
dicts the findings by Evstigneeva et al. (2006) who derived old 
ages for basically all UCDs in Virgo. 

It is clear from the discussion above that to decide whether 
dynamical masses of UCDs are really consistent with a pure 
stellar population is not straightforward. On the one hand, more 
accurate ages, metallicities and a-abundances are needed to 
better constrain the right choice of the stellar population mod- 
els (i.e. to constrain the IMF). On the other hand, there still is 
the freedom of adding an extended dark matter halo when mod- 
elling the UCDs for deriving their masses. In this case, also 
Bruzual & Chariot models with a Chabrier IMF could be the 
'right' choice. 

The most important result that comes out from our compar- 
ison between dynamical and stellar population model masses 
is that all UCDs and the nucleus of FCC 303 can, in principle, 
be explained by a pure stellar population. Dark matter is not 
mandatory for any of the objects within 1-3 half-mass radii. As 
mentioned before, it can not be ruled out, however, that UCDs 
are dominated by dark matter at large radii where high signal- 
to-noise spectra hardly can be obtained due to the very low 
surface brightness. 

5. Summary and Conclusions 

In this paper we have presented the spectroscopic analysis 
of four ultra-compact dwarf galaxies (UCDs) and one nu- 
cleus of a dwarf elliptical. These UCDs are among the five 



brightest in the central region of the Fornax cluster (Hilker 
et al. 119991 Drinkwater et al. 120031) . The analysis is based on 
high resolution spectra obtained with the instrument UVES 
mounted in the VLT at El Paranal, Chile. Velocity dispersions 
were derived from selected wavelength regions using the 
direct-fitting method by van der Marel (van der Marel & Franx 
119931 ). To derive the masses of the UCDs a new modelling 
program has been developed that allows to choose between 
different representations of the surface brightness profile (i.e. 
Nuker, Sersic or King laws) and corrects the observed velocity 
dispersions for observational parameters (i.e. seeing, slit size). 
The light profile parameters of the UCDs and the nucleus of 
FCC 303 were obtained from Hubble Space Telescope imaging 
(see the accompanying paper by Evstigneeva et al. [2006). The 
derived dynamical masses then have been compared to those 
expected from stellar population models. The main results of 
our analysis are as follows: 

(i) The observed velocity dispersions of the Fornax UCDs 
range between 22 and 30 kms -1 , slightly smaller than those of 
Virgo UCDs (Evstigneeva et al. 2006) and DGTOs (Hasegan 
et al. 2005). The observed velocity dispersions of the Na dou- 
blet regions are systematically higher than those of the other 
regions, probably reflecting the contribution of stars from dif- 
ferent evolutionary stages. 

(ii) The generalized King models for the one-component 
fits and King+Sersic models for the two-component fits were 
adopted to determine the masses of the UCDs, which range 
between 1.8 and 9.5 x 1O 7 M . The corresponding central and 
global projected velocity dispersions are in the range 29 to 
41 and 19 to 27 kms -1 , respectively. The masses as derived 
from Nuker and Sersic models depend on the adopted trun- 
cation radius and the slope of the outer power law (Nuker). 
Truncation radii that are too small lead to underestimated 
masses (e.g. Sersic model of UCD4) and flat outer power laws 
to a steeply increasing cumulative mass profile (e.g. Nuker 
model of UCD4). The central velocity dispersion mainly de- 
pends on the shape of the light profile within the innermost 
parsecs. If the cuspy behaviour of some UCDs is true, the cen- 
tral velocity dispersion would be up to 10 kms -1 higher than 
derived from generalized King models. 

(iii) The mass-to-light ratios MjLy of the Fornax UCDs 
range between 3 and 5. These values are compatible with 
predictions from stellar population models, especially those 
from Maraston (120051 ) assuming a Kroupa or Salpeter IMF. 
A Salpeter IMF, however, would imply that about half of the 
Fornax UCDs would be dominated by intermediate age pop- 
ulations of 4 to 8 Gyr. The high Balmer HyS indices of three 
UCDs in Fornax is consistent with this interpretation (Mieske 
et al. 120061 ). Also the M/L v values of all Virgo UCDs can be 
explained by stellar population models (see Evstigneeva et al. 
2006). No dark matter component is needed for UCDs whithin 
1-3 half-mass radii. One can, however, not rule out that UCDs 
might be dominated by dark matter in the outer parts. 

(iv) The mass-to-light ratio MjLy of the nucleus of 
FCC 303 is 2.5, entirely compatible with a pure stellar popula- 
tion and similar to those of other dwarf elliptical nuclei (Geha 
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et al. 120021 . 

We conclude that ultra-compact dwarf galaxies most prob- 
ably are the result of cluster formation processes (i.e. like 
large globular clusters, assembled star cluster complexes, nu- 
clear star clusters) rather than being genuine cosmological sub- 
structures themselves (i.e. compact galaxies formed in small, 
compact dark matter halos). 

Massive star clusters (and star cluster complexes) nowa- 
days form in interacting galaxies like the Antennae (e.g. 
Whitmore et al. U999ft . It has been shown that the most massive 
and compact young massive clusters (YMCs) and star cluster 
complexes can survive their 'infant mortality' (e.g. Fellhauer 
& Rroupa 2002 Bastian et al. 1200 6a). The most massive sur- 
viving YMCs have masses comparable to those of the UCDs 
presented in this paper (NGC7252:W3, Maraston et al. 120041 
NGC7252:W30, NGC1316:G1 14, Bastian et al. I2006bt . It is 
reasonable to assume that similar massive clusters have formed 
in merger events of the early epochs of galaxy cluster forma- 
tion. Some of them might have survived the disruption and 
merging of galaxies during the galaxy cluster formation pro- 
cess as free-floating 'UCDs'. 

But also the threshing of nucleated dwarf galaxies can not 
be ruled out definitely as formation process of UCDs. Nuclear 
star clusters as massive as UCDs were able to form either by 
the merging of globular clusters through tidal friction in early- 
type dwarf galaxies (e.g. Oh & Lin 120001 Lotz et al. 120011) or 
through successive central starbursts in late-type bulge-less spi- 
rals. During the dynamical evolution of galaxy clusters many 
of these low mass nucleated galaxies may be disrupted (e.g. 
Moore et al. 1999). Their debris then is distributed in the intra- 
cluster medium, building up the extended cD halos of central 
galaxies (e.g. Hilker et al. [1999). Their globular clusters and 
nuclei are expected to survive as intra-cluster stellar popula- 
tions, among them the present-day dark matter-free 'UCDs'. 

In order to differentiate between the formation of UCDs via 
super-star cluster formation in starbursts and the formation via 
tidal disruption of nucleated (dwarf) galaxies detailed analy- 
sis of the chemical composition and age structure of the stel- 
lar population in UCDs might help. A super-solar [a/Fe] abun- 
dance (as suggested for most of the Virgo UCDs, Evstigneeva 
et al. 2006) would point to a rapid, short-timescale formation 
process which is typical for old globular clusters. A galaxy nu- 
cleus would only show a super-solar [a/Fe] value if its star for- 
mation was truncated soon after its initial burst or if it would be 
composed of merged GCs which themselves have super-solar 
a abundances. However, present-day nuclei of dwarf ellipticals 
seemed to have formed in a rather continuous way since their 
a-abundances are enriched to solar values (Geha et al. 20Q3). It 
has to be seen whether the a-abundances of the Fornax UCDs 
are the same as for the Virgo UCDs. Intermediate age popula- 
tions in UCDs would indicate that some Gyr ago gas still was 
present in UCDs, constraining the threshing timescale and/or 
the time of the last major star formation event. Note that there 
are some hints that Fornax UCDs might be different to Virgo 
UCDs (e.g. Mieske et al. 120061 1. Maybe also the environment 
plays a role in their formation processes. 



Independent of which scenario is right it is interesting to 
note that star clusters more massive than about 5 x 1O 6 M seem 
to behave different than less massive 'ordinary' GCs. Not only 
the scaling and fundamental plane relations of massive clus- 
ters start to deviate from those of normal GCs (see results and 
discussions in Kissler-Patig et al. 120061 and Evstigneeva et al. 
2006), but also their stellar populations are getting more com- 
plex. Massive clusters do not possess a single stellar popula- 
tion, but are composed of sub-populations if different metallic- 
ities and/or ages. The best example is the most massive globular 
cluster in our Milky Way, a> Centauri, whose stellar populations 
revealed a broad metallicity distribution as well as an extended 
star formation history (e.g. Hilker & Richtler 2000 Hilker et 
al. 120041 1. Also Gl, one of the most massive GCs of M31 is 
suspected to contain multiple stellar populations (e.g. Meylan 
et al. 2001). It would be interesting to know whether this is true 
for all massive clusters, including UCDs. 

Although our study has shown that UCDs most proba- 
bly can be explained by pure stellar populations, more con- 
straints from observations are needed in order to get further 
insights into their nature. High signal-to-noise spectra for the 
Fornax UCDs would be useful to derive more accurate chem- 
ical abundances and age estimates from appropriate line in- 
dices. This should clarify whether they are enhanced in a el- 
ements and whether they posses intermediate age populations. 
Also high resolution spectroscopy (to obtain internal velocity 
dispersions) of more UCDs and galactic nuclei is needed to 
fill the parameter space (fundamental plane) in the interface 
between globular clusters and dwarf galaxies. The most ex- 
tended UCDs might be targets for spatially resolved high res- 
olution spectroscopy. This would show whether our modelled 
velocity dispersion profiles are real or whether one would need 
more sophisticated modelling to comply with possible orbital 
anisotropies and/or violation of the virial equilibrium through 
tidal effects. Spatially resolved velocity dispersion profiles also 
will allow to constrain the shape and mass of extended dark 
matter halos around UCDs if they existed. 
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